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In this lecture, an overview on interactions between atoms and surfaces is 
given that are mediated by the electromagnetic field. The emphasis is on dis- 
persion (or van der Waals) forces and transitions induced by thermal fluctu- 
ations in the near field of a surface. Applications with cold atoms held in 
microscopic traps near surfaces, as outlined in the lecture by J. Schmiedmayer, 
provide the experimental background for this theory lecture. 

We do not cover the chemical physics aspects that are relevant at distances 
comparable to the atomic scale. Phenomena like adsorption and diffusion on 
surfaces, or atomic beam diffraction from crystalline surfaces, are left to the 
lectures by J. Frenken and H. Lezec. 

43 : 

Outline 

> ■ 

First part: forces 

o : 

• Theory of the Van der Waals interaction. Electric dipole coupling, force 
from electromagnetic field fluctuations. 

o: 

• Characterization of the field near a surface. Cross spectral density, 
• ; fluctuation-dissipation theorem, example: planar surface. 

>>: 

43 . • Discussion of the Van der Waals-Casimir force. Asymptotic calculation 

of the near field spectral density. 



Second part: transitions 



Quantum states relevant for atoms 

Fermi's Golden Rule, connection to the cross spectral density 
Rates for different processes, scaling laws. 
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1 Van der Waals Forces 



1.1 Conventional viewpoint 

The Van der Waals force derives from the energy shift an atomic ground state 
experiences in front of a surface. This shift is computed in second order per- 
turbation theory, coupling the atom to the mode continuum of the electromag- 
netic field. The following section (in small characters like these) is taken from 
the lecture notes "Theoretical Quantum Optics I", held at Universitat Pots- 
dam (Germany) winter semester 2003/04. The complete notes are on line at 
www . quantum . physik . uni-pot sdam . de. 

1.1.1 Single-mode vacuum shift 

Let us start with a single mode of the electromagnetic field. We have seen previously 
that the ground state \g; 0) of the atom+field system is not affected by the interaction. 
This is actually only true in the rotating wave approximation. When we include the 
nonresonant terms 1 , 

H m = %g((7 + a) + er_a) 
they do affect the ground state, since 

H m \g;0) = hg\e; 1). 

In second order perturbation theory, we therefore get the following energy shift 

|(e;l|g nr |g;0)| 2 _ ft 2 g 2 
90 E el -E g0 Hcoa + u) W 

This shift is very small compared to those due to the resonant interaction because here 
the "detuning" uja + w is large (the transition can only happen "virtually" because it 
violates energy conservation). 

Nevertheless, we learn that the atom-field interaction does displace the absolute 
ground state of the system, even for a single mode. The Lamb shift is the generalisa- 
tion of this result for the full, multi-mode electromagnetic field. 



1.1.2 Lamb shift, a first glimpse 

The calculation of the Lamb shift is a basic example of the different infinities that occur 
in quantum electrodynamics. We only give a first idea of the calculation, without 

l g is a coupling constant with dimensions frequency, a and a* are the annihilation and 
creation operators for a photon in the field mode, a and cr* are the corresponding atomic 
ladder operators. 
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going into the many technical details (renormalisation, subtraction of classical energy 
shifts like Coulomb and polarisation self-interaction etc.). 

For the multi-mode field, any one-photon state \e; 1^) is coupled via the nonres- 
onant interaction to the ground state \g;0), giving a coupling matrix element hg^. 
Summing over all these modes, we find the energy shift 



—J2 



£f|d-s k/ J 2 |A M (x) 



(2) 



/ d3 S £0 nt + M ( |d|2 - |d - fe|2 ) |/k(x)p 



h (2vr) 3 7 2e V r (u; J 4 + kc) 

We have introduced the mode functions /k^(x) (equal to exp(ik • x) for a plane wave 
expansion) and, assuming that | /k (x) | 2 is independent of /U, performed the summation 
over the polarization unit vectors (as already seen a few times in the exercises). Good 
luck, the volume V cancels with the factor V -1 / 2 from the 'field per photon' E k . The 
integration over the angles of k gives: 

J sm6&6&v (|d| 2 - |d • k| 2 ) = y|d| 2 . (3) 

We are finally left with the radial integral 

8^ Idl 2 /" A k 3 dk 



6 E g0 = — 



d f k dk < \ 

— / T~TT. wr ong (4) 

e Jo ua/c + k 



6(2tt) 2 

This integral is obviously divergent in the limit k — ► oo, which is called an "ultraviolet 
catastrophe". The integral is only finite if we introduce a cutoff wavenumber A. What 
order of magnitude can we reasonably give to A? Remember that our whole theory 
is based on the "long-wavelength approximation" (photon wavelength 2ir/k much 
larger than atom size ~ do). For wave vectors k S> 1/ao we therefore cannot trust our 
interaction Hamiltonian any more. Choosing the cutoff A = 1/ao, w e get a shift of the 
order of d 2 A 3 /e = d 2 /ea^ ~ e 2 /eoao- This is a completely wrong result because it is of 
the order of the atomic binding energy! 

To get the right result, one has to take into account self -energy terms that we did 
not explicitly write down in the Hamiltonian. These are infinite, too, but if they are 
evaluated with a similar cutoff and subtracted, one gets something which does not too 
badly diverge with the cutoff A: 

5E 90 « in (^V) ( c ° rrect ) < 5 ) 

3ir 2 Eqc 6 \uja/cJ 

Here, the correct value of the cutoff is given by . . . the Compton wavevector of the 
electron, A = mc/h (up to this wavevector is our nonrelativistic theory valid). 
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It is beyond the scope of this lecture to give a more detailed account of the Lamb 
shift calculation. 



1.1.3 Van der Waals and Casimir-Polder forces 

The previous calculation, as sick as it is, does give a result for the modification of 
the Lamb shift due to boundary conditions imposed on the electromagnetic field by 
macroscopic bodies. Let us focus on the simple case that our two-level atom is placed 
at a distance z from a perfectly reflecting mirror. All we have to change in the calcula- 
tion are the mode functions, and we shall use 

/k M (x) = V2sin(k z z) exp[i(k x x + k y y)}. 

(We are actually cheating with the polarisation vectors to simplify things.) The factor 
\/2 comes from the fact that we use a unitary transformation 2 of the plane waves 
e ±ifc z 2 ^y e have |- change the angular integration {3)l because now the z-axis plays a 
privileged role. We can expand the dipole moment in cartesian components and get 
(9 is now the angle between k and the z-axis): 

"i2 f 2 

d<p |d • k| = / dip d x sin 9 cos p + d y sin 9 sin p + d z cos 6 



2tt (\d 2 x sin 2 9 + \d 2 y sin 2 9 + d\ cos 2 



Note that the integrals over the mixed terms involve cos p, sin <p, or cos <p sin <p whose 
average over one period vanishes. In terms of the components d\\ = d z and d± of the 
dipole parallel and perpendicular to the z-axis, we thus get 

1 fk 3 dk sin6>d6> 

oE g0 = 



2(2ir) 2 e J u; a /c + k 
x (£ sin 2 9 + d\(l - I sin 2 #)) (1 - cos(2A;z cos 9)) (6) 

The '1' in the last bracket is independent of the distance z to the mirror — it will 
give the Lamb shift for an atom in free space. The second term with the cos gives the 
modification of the Lamb shift we are interested in. The 9 integral gives (make the 
substitution u = cos 9 and differentiate with respect to 2kz to get the integrals with 
the additional factor cos 2 9) 



sin 9 d9 cos(2kz cos 9) 



sin 9 d9 cos 2 9 cos(2kz cos 9) 



2 sin 2kz 
2kz 

2s'm2kz 4cos2kz 4sin2kz 



o y ' 2kz (2kz) 2 (2kz) 

f w . . . 2 . „. 4cos2/cz 4sin2fcz 

/ sm9d9 sin 9cos(2kzcos9) = — — — h 

Jo (2kz) z 



(2kz) 



2 Alternative argument: normalise the volume integral of the squared mode function to V. 
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Apparently, we get an expression with different powers of 1/z. 

Let us focus first on the short-distance limit and consider only the "most diverg- 
ing" term ~ 1/z 3 . It involves the integral 



Hz) 



00 dk sin 2kz 
kA + k 



where we have put kA = oja/c. We use a trick in the complex plane to perform the 
integration: write sin 2kz as the imaginary part of e 2lkz and deform the integration 
path in the complex plane from the positive real axis to the positive imaginary axis. 
This works because the exponential e 2lkz vanishes at infinity, and the integrand has no 
pole in the upper left quadrant. Using the integration variable k = in on the imaginary 
axis, we find 

fQ o dke 2ikz f°°dKe~ 2KZ 
h\z) = Im / —j- = Imi 

= k A 



kA + k Jo kA + ik 

°°dKe- 2KZ 



Jo k\ + K 2 ( } 

We approximate this integral by noticing that it is the product of two decaying 
functions. The exponential decays on a scale given by 1/z, while the "Lorentzian" 
l/{k\ + k 2 ) decays on the scale kA- Now look at which scale is larger. The short dis- 
tance regime is defined by the limit z <C the distance is much smaller than the 
transition wavelength. Then the exponential decays much more slowly, and we can 
replace it by unity. This gives 



kAZ <C 1 : H z ) = kA [ , o ^ o = arctan(K,/kA 

Jo ki + K 



00 7f 

~ 2 



The other inverse powers can be handled in a similar way and give sub-dominant 
contributions. We thus find the so-called van der Waals energy 

E vd w = --. 11 \ - {z < Xa, wrong). (8) 

47T£o 0Z 6 

This is quite close to the correct result. The correct calculation with the proper hand- 
ling of the polarisation vectors gives 

1 df + K 

#vdw = ~ 47r£o 8z 3 — (z -^Xa, correct). (9) 

which differs only in the sign of the d\ term. For realistic atoms, you get contribu- 
tions from all excited states that are connected to the ground state (and from all their 
Zeeman sublevels). 
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In the long-distance limit k A z 3> 1 (distance larger than transition wavelength), 
we can also evaluate the integral Q. Now the exponential has the fastest decay, and 
we can replace the lorentzian by a constant: 
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k A z > 1 : h(z) = — I dree 



— 2kz 



k A Jo 2k A z 

and hence a result much smaller than in the short distance limit (because now k A z ^> 
1). This occurs in fact for all terms we got. The final result is the famous Casimir- 
Polder energy: 

Scp= -idb? {z>>Xa) - 

where d 2 = cZ| + d\ (no distinction between the dipole orientations). 
1.1.4 Summary 

An atom in its ground state close to a planar surface feels an attractive po- 
tential that scales like 1/z 3 (van der Waals potential). This applies to distances 
much shorter than the wavelength of the strongest electric dipole transition 
starting from the ground state. The coefficient depends on the squared dipole 
moment d 2 of the atom, more precisely its expectation value in the ground 
state. At large distances z 3> \ A/ a power law l/z 4 is found (Casimir-Polder 
potential). The coefficient depends on d 2 /uj A , and this quantity is proportional 
to the atom's static polarizability a(0). More details on this calculation can be 
found, for example, in the textbook by Craig and Thirunamachandran |1] and 
the 1990 Les Houches lecture by Haroche [2J. 

For later reference, we quote here the result for the polarizability of the 
atomic ground state. This formula can be found by evaluating to lowest order 
in an applied electric field (frequency cu) the average dipole induced moment. 

/ ^ uj eg {g\d i \e){e\d j \g) 
aiM ~h\ c^-^-Oic (10) 

Here, the Bohr frequency between the ground state g and some excited state 
e is uj eg = (E e — E g )/h. The infinitesimal imaginary term in the denominator 
ensures causality: as a physical response function, the polarizability ^-(o;) 
can only have poles in the lower half of the complex cu-plane. Note that a 
number of different conventions for the units of the polarizability occur in the 
literature. In our notation, a^- (a>)/eo has the units of volume, and the average 
induced dipole moment is given by {di)e~ luJt + c.c. = otij{u)Ej{r) e~ ltJt + c.c. 
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(Here and in the following, we use Einstein's summation convention and drop 
the sum •.) 

For the static polarizability, we get 

, m _ 2(p|di|e)(e|dj|ff) 

If the ground state has spherical symmetry, this tensor is actually proportional 
to Sij. (Note that in the hydrogen atom, this requires the summation over the 
three p x ,p y , p z excited states.) 



1.2 Role of field fluctuations 

We now present an alternative viewpoint on the van der Waals force. Its ad- 
vantage is that the role of field fluctuations is made explicit. Both vacuum and 
thermal fluctuations contribute at finite temperature. We develop a formalism 
that treats both in a unified way. At zero temperature, the standard van der 
Waals force is recovered. A thermal correction appears in the near field of a 
'hot' surface. 

The following section is adapted from the paper "Radiation forces on small 
particles in thermal near fields" by C. Henkel, K. Joulain, J.-P. Mulet, and J.-J. 
Greffet [J. Opt. A: Pure Appl. Opt. 4 (2002) S109]. 



1.2.1 Radiation force 

We consider a model system made from a dielectric half-space and a small (sub- 
wavelength) "test particle" in vacuum above. Both objects are possibly at finite tem- 
perature and radiate electromagnetic fields whose sources are thermal current fluctu- 
ations. Using the fluctuation electrodynamics framework of Lifshitz |3|, the spectral 
density of these currents is related to the imaginary part of the dielectric function (for 
the half-space) or of the polarizability (for the particle). 

We are interested in the force exerted by the radiation on the test particle. This 
force is actually both a time-average over the rapid fluctuations of the thermal field 
and an average over a statistical ensemble for the field. For particles with a subwave- 
length size, the force is given by 1 4 1 

F(r)= J2 (di(t)VEi(r,t)) (11) 

i=x,y,z 

where d(t) is the dipole moment of the particle and E(r, t) the total electric field at the 
particle's position. The latter contains the fields generated by the thermal substrate 
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and by the particle's dipole moment, as well as the blackbody field in the vacuum 
half-space. The expression iTTl combines both the Coulomb force (corresponding to 
the potential energy — d • E) and the Lorentz force (involving the magnetic field), see 
l4l . As a simple derivation, we note that the force operator can be defined in terms the 
Heisenberg equation for the momentum operator 

P = p=^[H,p] = -VH int (12) 

where we have used the fact that the atomic Hamiltonian is independent of the atomic 
center-of-mass position. Using the electric dipole interaction, one gets fill . 

In our case, both the dipole moment and the electric field are fluctuating quantities. 
The force may therefore be written as a sum of two terms 

F(r) = (4 nd \t)VEf\r,t)) + (df ) (t)V J ef nd) (r,t)) (13) 

where the first describes the (spontaneous and thermal) fluctuations of the field that 
correlate with the corresponding induced dipole, while the second involves dipole 
fluctuations and the field they induce. There is no term involving fluctuations of both 
the dipole and the field: these are not correlated, since they originate from different 
physical systems. More precisely, to lowest order in perturbation theory, we can take 
the fluctuations in (I3j to be those of the non-coupled atom+field system, these being 
obviously uncorrected. 

1.2.2 Induced fluctuations 

The dipole induced by the field fluctuations is given by the particle's polarizability: 

d( ind »H=a( W )E(r; U ) (14) 

where we switched to frequency Fourier transforms and assumed an isotropic polar- 
izability. Note also that E(a>; r) is the total electric field at the dipole's location. If we 
work at lowest order in the polarizability, we can identify this field with the fluctuat- 
ing field E( fl ) (r; uj) and ignore the scattering of this field by the particle. Otherwise, we 
could work with a "dressed" polarizability. Note however that the fluctuating field is 
not simply a collection of plane waves, but takes into account their scattering from the 
substrate. 

The field induced by the dipole fluctuations is given by the electromagnetic Green 
tensor: 

E( ind )(x;w) =G(x,t;lj) ■ d(o/), (15) 

where the variables x and r represent the observation and source point, respectively. 
To lowest order, we can identify the total dipole in llT5l with its fluctuating component, 
as we did before for the field. 
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Finally, we express the operator product in terms of frequency Fourier transforms 

as 

(A(t)B(t)) = J ^e^-«'»(A\u)B(J)) (16) 

where the frequency integrals run over both positive and negative frequencies. With- 
out loss of generality we have written the frequency integral for B(t) in terms of the 
conjugate operator £?t[cj]. This is possible because in < fl3b . both A(t) and B(t) are real 
(hermitian) operators. 

1.2.3 Fluctuation spectra 

Dipole. The dipole fluctuations of the particle are characterized by a spectral den- 
sity given by the fluctuation-dissipation theorem. This assumes the particle to be at 
thermal equilibrium (temperature T2/&.B) and gives 0|6| 

<4(«KV)> = 2vr5( W - "Oftj ^/T^ W") (17) 

Note that this expression applies to both positive and negative frequencies: for uj > 0, 
we have ( e hu} / T 2 — 1) _1 = n{uj/T2), the average excitation number for Bose-Einstein 
statistics (vanishing at T = 0). At negative frequencies, uj < 0, we have (e hu3 / T2 — 
l)- 1 = -(n(M/T 2 ) + l),andthe minus sign is compensated by the symmetry relation 
Im OLij(uS) = — Im a*ij [—oj] so that the total expression is positive again (as it should). In 
particular, at T 2 = 0, the dipole "fluctuates only at negative frequencies" which reflects 
the existence of the ground state (fluctuations can only connect to higher-lying states). 
We identify these fluctuations with spontaneous fluctuations. Thermal fluctuations 
are proportional to n and are the same for positive and negative frequencies (as they 
should for classical noise spectra). 

Field. For the field fluctuations, the fluctuation-dissipation theorem reads (field at 
temperature T x ) EJ|6| 

(^(rsu^rV)) = 2vr^-^) eWTi _ 1 ImG tJ (r,rV) (18) 

where the imaginary part of the Green tensor appears. This result is valid at equi- 
librium (field and source at the same temperature). For a generalization to non- 
equilibrium, see below. 

1.2.4 Resulting force 

Collecting the different contributions to the force, we show that the "vacuum" fluctu- 
ations of both particle and substrate conspire to give the standard van der Waals- 
Casimir-Polder-London force. The thermal fluctuations give a correction that we 
evaluate in the limit of a particle at zero temperature. 
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Starting from the general expression il3l . the force due to dipole fluctuations is 
(using the notation cL; = &u/2it) 

(df\t)VEt d \r,t)} = 2h f ( jj£_ 1 (Ima(u;)) ViG«(r,r ;tJ ) (19) 

where Vi is shorthand for the gradient with respect to the first position variable of the 
Green tensor. 

The force due to field fluctuations is 

(df d) (i)v4 fl) (r^)) vac = 2h J&> J e ^_ 1 «*H (ImV 2 G ii (r,r; W )) (20) 

We now argue that the Green tensor satisfies the symmetry relation 

ViG^r, r; u) = V 2 G u (r, r; u). (21) 

Very generally the reciprocity theorem implies that the Green tensor is symmetric 
under combined exchange of the position arguments and vector indices |5|. For the 
particular case of the half-space, we have checked this explicitly. In this case, we also 
find that only the part involving the reflection from the interface remains under the 
derivative. This is because the (imaginary part of) the free space Green tensor depends 
only on (r — r') 2 , due to isotropy and has vanishing first derivative at r = r'. 

Using these results and specializing to the equilibrium case T\ = T 2 = T, we can 
combine Jl9l > . ODl to give the van der Waals force 

f p u/T _ l a(w) ViGa(r,T;u). (22) 

You will find frequently in the literature an expression where the u-integral is trans- 
formed in the complex u;-plane. One closes the integration path with a semicir- 
cle in the upper half plane. This circle does not contribute to the integral because 
a(uj) and Gij(u), being physical response functions, decay to zero for \u>\ — > oo. 
In addition, they have no singularities in the upper half plane. The integral along 
the closed contour is thus determined by the poles of the Bose-Einstein factor at 
uj = i£ n = 27rinT/h, n = 0,1,2... ('Matsubara frequencies'). The pole n = con- 
tributes with only one half of its residue because it lies on the original contour. One 
thus gets: 

F vdW (r) = 2T£'a(i£ n ) ViGii(r,r;i£ n ). (23) 

n 

(The prime reminds that n = is counted with a factor ^.) We have dropped the Tm' 
because the response functions are real at purely imaginary frequencies. 
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Zero temperature. In the limit of a particle at zero temperature, the Bose-Einstein 
distribution becomes a step function, 

lim — -,= = H(-uj), (24) 

showing that the integral runs only over negative frequencies. We thus find 

rod 

Fvac(r) = 2Mm / &u a(u) ViG^(r, r; u), (25) 
Jo 

where we have used the symmetry properties of a and G%j . This result resembles up to 
a factor of 2 and a — Vi the expression for the van der Waals-Casimir-Polder-London 
energy given in |6|. The gradient occurs because we calculate the force, of course. To 
understand the factor 2, we note that when the force is calculated from the energy, 
both position dependences in Gu(r,r' = r;uj) get differentiated. Close to a planar 
surface, both give the same contribution, since in the Weyl plane-wave expansion (see 
I1.4.H , the reflected part of the Green tensor is proportional to e n ^ z+z '^ where 7 is the 
normal component of the wave vector. In the planar geometry we consider here, the 
force is oriented perpendicular to the interface, of course. 

We finally note that it is a priori not clear that the light force derives from a potential 
because we are dealing with a dissipative system. At finite temperature and in equi- 
librium, the force is related to the energy change when the particle is displaced. This 
creates electromagnetic fields that do work on the polarization fluctuations, which is 
partly dissipated in the system. The corresponding potential energy is given by the 
free energy of the system. More details can be found in |7| and references therein. 
Here, at zero temperature and at lowest order in the polarizability, we do find a po- 
tential since F can be written as a gradient. This need not be true at higher orders, 
however, because the combined state of atom and field develops correlations (the den- 
sity matrix contains a term exp(d • E(r)/T) involving the interaction Hamiltonian). 



1.3 Fluctuation-dissipation theorem 

We outline here a proof of the fluctuation-dissipation theorem for the electro- 
magnetic field close to some arbitary macroscopic object. The only require- 
ment is that the object interacts in a linear way with the field and that ther- 
mal equilibrium prevails. This text is adapted from the review paper "Electro- 
magnetic field fluctuations on the nanometer scale" by C. Henkel, submitted 
for publication in the Handbook of Theoretical and Computational Nanotechnology 
(2004). 

The fluctuations of the electromagnetic fields can be characterized by statistical 
electrodynamics, which results from the application of equilibrium thermodynamics 
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and quantum theory to the macroscopic Maxwell equations |8|. We introduce in this 
Section the basic definitions for the field fluctuation spectra and review how they can 
be calculated. A key result in this context is the fluctuation dissipation theorem de- 
rived by Callen and Welton |9| whose proof is sketched here. We conclude with some 
remarks on how to handle non-equilibrium situations. 



1.3.1 Basic definition 

We shall assume that at thermodynamic equilibrium at temperature T, the field and 
the solid medium can be described by a Gibbs ensemble: each state with energy E is 
weighted with the Boltzmann factor ex.p(—E/ksT) = e~@ E . In the quantized version 
of the theory, these states are stationary states: they are eigenstates with energy E 
of the corresponding Hamilton operator H. We shall work in the Heisenberg picture 
where the field variables (called observables) evolve in time, while the state of the 
system is fixed. The Gibbs ensemble is then characterized by the density operator 

exp(-H/k B T) 
Pe q = ~ g— — ■ (26) 

Tr exp(— H/kbT) 

This is an operator on the Hilbert space of the system that can be represented by a 
(infinite-dimensional) matrix, the density matrix. Tr exp(— H/ksT) is the trace of the 
operator-valued Boltzmann factor, it is also called the partition function |10|. 

With respect to this equilibrium ensemble, we can define average values for the 
observables of interest. This average combines the quantum expectation value in a 
given stationary state with the corresponding statistical ensemble weights. The aver- 
age electric field, for example is given by 

(E(x, t)) = Tr [E(x, t)p cq ] = Tr [peqE(x, t)} , (27) 

where we have made use of the cyclic permutation under the trace. The time- 
dependence of the field is generated by the Hamilton operator, so that we have 

(E(x, t)) = Tr [exp(ifft/n)E(x) exp(-iHt/h)p cq \ = Tr [E(x)p eq ] , (28) 

where E(x) is the electric field operator at time zero and we have used that the Gibbs 
density operator i2"6l is invariant under time evolution. The average can now be com- 
puted in the Schrodinger picture and is found to vanish at equilibrium. In the classical 
theory, this is because the phase of the field is uniformly distributed. In the quantum 
theory, the stationary states for each mode of the field (labelled by k) are eigenstates of 
the photon number operator a) K a K . The field observable is a linear combination of the 
annihilation and creation operators a}, and a K that lower or raise the photon number: 
their quantum expectation values thus vanish in a photon number eigenstate. More 
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details can be found in Section !! .3.21 and in the textbooks by Loudon 1 11 1 and Mandel 
and Wolf 113 . 

The relevant information about the field fluctuations is thus encoded in the corre- 
lation function 

<E(x,t)E(x',t')> =Tr [E(x,t)E(x , ,t')/5 cq ] = (E(x, 0)E(x', t' - i)>. (29) 

In the second step, we have shifted the time arguments using the fact that the time 
evolution commutes with the equilibrium density operator ll26l . As expected from 
stationarity this correlation function only depends on the time difference r = t' — t. 
In the limit |t| — * oo, one expects the fields E(x, t) and E(x', t') to decorrelate and the 
correlation function {29jl to vanish. The time scale on which this happens gives the 
correlation or coherence time of the field. 

The spectrum of the field fluctuations can be defined by the Fourier expansion of 
the correlation function 3 

+oo 

(E(x, t)E(x', t')) = J ^e-W-QSfr x'; u). (30) 

— oo 

This relation is also known as the Wiener-Khintchine theorem |12|. The spectrum is 
actually a tensor, the cross-spectral density tensor. For x = x', its components specify 
the local polarization state of the field. As a function of x— x', it characterizes the field's 
spatial coherence, i.e., the contrast of interference fringes in a double-slit experiment 
with slits placed by x and x', see, e.g., the textbooks by Mandel and Wolf |12| and by 
Goodman 1131 . As a function of the frequency lo, the spectrum specifies the strength 
of the field fluctuations. This can be seen from the equivalent relation for the Fourier 
transforms of the field 

(E(x, -w)E(x', u')) = 2k5(uj - lo') £(pc, x'; lo') = ([E(x, w)] t E(x / , u>')). (31) 

In the last step, we have made use of the reality of the electric field that leads to iden- 
tity [E(x, Lo)y = E(x, —lo) for real lo. This way of writing also shows that the diagonal 
elements of £ (x, x; lo) are positive. More generally, < u* • £(x, x; lo) • u for any com- 
plex vector u. We note that the Fourier transforms of the fields strictly speaking do not 
exist as ordinary functions. The Fourier calculus nevertheless applies symbolically for 
the corresponding operator-valued distributions. 

The fluctuation spectrum of the electric field plays a key role for spontaneous and 
stimulated decay on the electric dipole transitions of an atom or molecule. The cor- 
responding spectrum for the magnetic field characterizes the perturbation the field 
exerts on an atomic magnetic moment or spin. This is discussed in detail in Section^ 

3 Boldface vectors placed near each other mean a tensor product, often written ®. Tensors 
are written in calligraphic letters like EorH. 
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1.3.2 Blackbody fluctuations 



Let us illustrate the correlation function introduced above with the example of the 
blackbody radiation field in free space. The electric field operator can in that case be 
expanded in plane wave modes | i~TJIT2l 



E(x,t)=£ 



2e V 



[a K (t)£exp(ik-x) +H.C.] , 



(32) 



where the mode label k = (k, e) combines the wave vector k and the polarization 
vector e _L k, and V is the quantization volume. We assume periodic boundary con- 
ditions so that the allowed wave vectors are discrete. The mode frequency is given 
by u(k) = ck = cV\i*. 'H.C.' denotes the hermitean conjugate operator so that the 
electric field is globally hermitean. In the absence of any sources, the annihilation and 
creation operators evolve according to 

a K (t) = a K exp[— iu(k)t], a K (t) = «I exp[ico(k)t]. (33) 

The Schrodinger operators a K and aj. satisfy the bosonic commutation relations 



= a K 4' - a i> « K = S K, «' 



(34) 



The equilibrium expectation value of products of the mode operators decorrelates for 
different modes because the density operator factorizes into a product over all modes. 
Therefore, the expectation value (aL a K ) vanishes for k / k'. 

For a given mode k, one gets from the Gibbs ensemble the Bose-Einstein occupa- 
tion number (with f3 = l/(k B T)) 

(a j K a K ) = n[u(k)] 



l + n[uj(k)} 



e f3hui(k) _ 

1 



1 



-I3huj(k) ' 



(35) 
(36) 



while the products a K a K and a^aj. have zero average. To prove Eq. 051, we evaluate 
the trace in the eigenbasis \n K ) of the photon number operator a\a K for the given 
mode. The energy eigenvalue of \n K ) is huj(k)(n K + i), and the summation over the 
Boltzmann weights gives 



(4<0 



oo 
n K =0 



J2 e -Phuj(k)(n K + \) 
ra K =0 







^ lo g(E e " 



\n=0 



|log(l-e^) 



0( 



1 



(37) 
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Using these informations, a straightforward calculation in the continuum limit, 

E k V J d 3 fc/(2vr) 3 , leads to 

f(x,x + r;w = Y 1-fr V / 7 38 



V [uir/cY (Lor/c) 6 J ) 



which only depends on the difference vector r = x' — x, as expected. This expression is 
very similar to the Green tensor in free space which is not a coincidence, but a special 
case of the fluctuation-dissipation theorem introduced in Section lT.3.31 

In the limit r — > 0, one can check that £ (x, x; uj) is proportional to the unit ten- 
sor and positive. Taking the trace, one gets, up to a factor eq/2, the spectrum of the 
electric energy density u e (uj). Summing the contributions of positive and negative 
frequencies, the electric energy density is given by the Planck formula 

oo 

/du^ 
g-Wbb.eC^), ( 39 ) 



«bb,e(w) = ^ ^ = 2ttDOS (w)— n (w + - . (40) 

In the last expression, we have made use of the free space dos(cj) = uj 2 /tt 2 c 3 . This 
result has an intuitive explanation: the electric energy density is the density of modes 
per unit volume, DOS(u;)da;, multiplied by one half of the average equilibrium en- 
ergy Uu{n{uj) + |) of a harmonic oscillator (the other half contributes to the magnetic 
energy). 

At zero temperature, the electric plus magnetic energy per mode is given by the 
ground state oscillator energy huj/2, and this leads to a divergent integral in the UV. 
From this divergent zero-point energy, one can extract a finite, measurable energy 
difference, called the Casimir energy, when the mode functions are changed by the 
boundary conditions imposed by material structures. See the textbook by Mostepa- 
nenko and Trunov |14| for more details. The energy measured by a photodetector 
does not diverge because it is proportional to average of the photon number operator 
(a}.a K ) = n(u){k)) 1TT1 fl2l . If the creation and annihilation operators are ordered in 
this way, the divergent contribution of the zero-point energy disappears. The black- 
body spectrum then shows an exponential decrease beyond the thermal wavelength 
A t h = 2nfic I (&bT) (the Wien displacement law), and the frequency integral 0^1 be- 
comes convergent. 
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1.3.3 Fluctuation-dissipation theorem 

The previous example suggests that there is a relation between the field fluctuation 
spectrum and the Green function. With the previously introduced notation, we have 

S(x,x';u) = e p^_ 1 Img(x,x», (41) 

where the Green tensor is defined as the electric field radiated by a monochromatic 
point dipole, 

£f ip) (x; u) = 0y(x, r; w)d», (42) 

and the fluctuation spectrum by Eq. J30b . Eq. d4lTt is actually true under more general 
conditions and is known as a fluctuation-dissipation (FD) theorem |9[. It holds for 
linear systems and their fluctuations around the thermal equilibrium state. The dis- 
sipation is encoded in the imaginary part of the response function that characterizes 
the linear response of the system to an external perturbation. The FD theorem is of the 
form given here provided the Green tensor satisfies the symmetry condition specified 
in Eq. ll2"l"t . 

The FD theorem will be our basic tool to compute field fluctuation spectra in the 
near field of nanostructures. We give in this Section an introduction for linear sys- 
tems and summarize a general proof in the context of statistical electrodynamics. A 
generalization to non-equilibrium situations is discussed in Section li.3.41 

Johnson-Nyquist noise in metals. Consider a resistance at temperature T. One 
observes a thermal fluctuation of the current through the resistance, called Johnson 
noise, whose variance in a given bandwidth Af = Alo/(2it) is given by the Nyquist 
formula 

<" >a " = ^W) ■ (43) 

We now show that this formula can be related to a fluctuation dissipation theorem 
for the current density j(x) of the resistance. Consider first the z -component of the 
current density and a small volume element AV = Az AA. From the current noise 
along the ^-direction, one then has 

2 Aoj ik B T Aco 4k B Ta(u) 

{SM )A - = ^WjAA* = ^^X^' (44) 

where a(oo) = Az/(R(u)AA) is the conductivity (the inverse of the specific resistance). 
In the low frequency range where the Nyquist formula is valid, the conductivity can 
be expressed via the dielectric function of the resistance, eoe(uj) = eo + ia/u. The 
result |@4) can thus be obtained by averaging the following relation over the volume 
element 

(Jj(x)^j(x / ))A^ = -^Ak B Tu> llm£ £(x;a;)5(x - x'). (45) 
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We have assumed that neighboring volume elements have uncorrelated current noise, 
hence the spatial delta function. We use here the convention that the current noise J45l 
is given by the integral of the noise spectrum J7(x, x'; uj) over intervals Auj/(2tt) cen- 
tered at positive and negative frequencies ±u. Since o;Ime(x;Lj) is an even function 
of uj, the noise spectrum is given by 

J(x,x?;w) = 2k B Tu 1 Ime e(x; w)<5(x - x'). (46) 

This result has already the structure of the FD theorem <4H . Since the dielectric func- 
tion gives the polarization induced by an electric field, a more natural formulation is 
in terms of the polarization noise spectrum (writing j(w) = — iwP(w)) 

V(x,x';uj) = — Hm£ e(x;wWx-x'). (47) 

UJ 

The temperature-dependent prefactor is the low-frequency limit of 2h/(e hu) ^ kBT — 1) 
occurring in Eq. ll4Tll . The spatial 5-function applies to a local dielectric response. The 
generalization to a nonlocal medium is immediate: the fluctuations are then correlated 
on some characteristic scale, typically the mean free path. 

We thus find that the strength of the thermal Johnson noise at a given frequency in 
an absorbing material is related to the amount of dissipation, as encoded in the imag- 
inary part Im e(u). This permits to characterize the thermal polarization and magneti- 
zation fluctuations that appear in the macroscopic Maxwell equations. The noise spec- 
trum of the polarization noise Pfl(x; uj) is given by the Johnson-Nyquist formula j47} . 
with the factor k&T replaced by huj/(e^ huJ — 1) to be valid at all frequencies. If the ma- 
terial is magnetic with a (local) permeability /i(x; uj), it contains magnetization fluctu- 
ations Mfl(x; uj) with a spectrum 

Mx,x ; .) = -^l m; ^ 5 (x-x'). (48) 

The fluctuating material polarization radiates an electromagnetic field that, in thermo- 
dynamic equilibrium, compensates for the loss of electromagnetic energy inside the 
material. Only in this way is it possible to enforce the equipartition law of thermal 
equilibrium, every degree of freedom (here the polarization field) carrying an energy 
UbT/2. Consider the balance of electromagnetic energy for a system without external 
polarization sources. On the one hand, the mechanical work performed per unit time 
can be written as 

ujIm(P\(x;uj) ■ E(x;u/)) 

= 2ttuj5(uj — uj 1 ) J d 3 x' Vij(x, x'; uj) Im £/y(x, x'; uj) 

= 2tt5(uj — ^')—0f— ^hri£o£(x; uj) ImTr Q(x, x; uj). (49) 
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We have taken into account that only the part of the field radiated by the polarization 
fluctuation is correlated with the latter and expressed that field in terms of the Green 
tensor {42jl. A similar result holds for the magnetic contribution. On the other hand, 
the average over the electric losses in the medium leads to the same expression 

a;Im[e e(x;Lij)](E^(x;u;) -E(x;u/)) 

= 2tt5(uj - u/)Im[e e(x;ij)] - ImTrS(x,x;u)), (50) 

using the fluctuation-dissipation theorem l l4lTl . As a result, the energy the polarization 
emits into the field rEq.d49b1 is exactly compensated for by the field energy lost by 
absorption, Eq. dSOb . This also implies that the average Poynting vector Re(E^ x H) 
vanishes, since there is no net energy transfer between medium and field. 

We shall see that in the quantized theory the polarization fluctuations of the mate- 
rial also contribute to the quantum fluctuations of the field. Otherwise, the dissipation 
present in the macroscopic Maxwell equations would force the field operators to de- 
cay to zero, including their commutators. These are preserved due to the quantum 
fluctuations of the material polarization. This picture suggests as well the existence of 
an FD theorem: the material loss that forces the fields to decay must be balanced by 
the fluctuations inside the material. 

Properties of quantum field fluctuations. Before giving a general proof of the 
FD theorem i4Ti , let us summarize some of the properties it implies for the equilibrium 
fluctuations of quantized fields. 

The FD theorem allows to compute the quantum and thermal fluctuations of the 
electromagnetic field once the Green tensor is known. This quantity can be computed 
by solving the macroscopic Maxwell equations with point-like dipole sources. As 
long as the medium responds linearly to the field, the classical version of the theory is 
sufficient, quantum and thermal fluctuations are handled self-consistently using the 
theorem. 

The noise spectrum of a quantized field is not symmetric. It is proportional to the 
Bose-Einstein occupation number n(uS) for positive frequencies and decays to zero for 
hu 3> k^T. At negative frequencies, one finds, given that ImQ(u) is an odd function, 
that the spectrum is proportional to l+n(|aj|). The zero-point fluctuations appear here. 
The asymmetric frequency spectrum of zero-point or vacuum fluctuations can be un- 
derstood qualitatively by noting that in the ground state, a system can only fluctuate 
via a virtual transition towards a state with higher energy. The corresponding Bohr 
frequencies are all positive. (That this leads to a spectral weight at negative u is related 
to our - conventional - choice of the exponential factor in the noise spectrum lOOl . - ) In 
the high-temperature limit or, equivalently, for classical systems, upward and down- 
ward transitions occur with equal probability, and the fluctuation spectrum is sym- 
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metric: n(u) ™(M) + 1 ~ k^T '/Tuo 3> 1. For intermediate temperatures, we show 
below that the principle of detailed balance is satisfied, up- and downward transition 
rates differing by a factor e@ huj . 

At positive frequencies, Eq. J30l> shows that the spectrum £ (x, x'; oj) picks out that 
part E( + )(x',£') of the electric field operator that evolves like e~ lu)t ' . By analogy to 
time-dependent wave functions in quantum mechanics, this part is called the positive 
frequency part of the field. In the mode expansion ll32l . it corresponds to the sum over 
the annihilation operators a K [see also Eq. ll33l 1. Similarly, only the negative frequency 
part E( - )(x, t) = [E( + )(x, of the field operator contributes in the first factor. It 
follows that in the vacuum state, the expectation value 



Tr 



E(-)(x,t)E( + )(x',t , )Pvacl =0 (51) 



vanishes since the annihilation operators, by definition, give zero when acting on 
the vacuum state. This operator order (annihilation operators acting first) is usually 
called 'normal' order. A typical example is the intensity measured by a photodetec- 
tor li TTllT2l . A nonzero vacuum expectation value occurs with the reverse operator 
order (creation operators acting first). This anti-normal order is picked out for nega- 
tive frequencies in the fluctuation spectrum £ (x, x'; u>). It gives nonzero results even 
in the vacuum state (at zero temperature), as we have seen in Eq. i36l and in the FD 
theorem d4lT> . 



Proof of the FD theorem with linear response theory. As mentioned previ- 
ously, the macroscopic Maxwell equations in an absorbing medium have to be sup- 
plemented by material fluctuations in order to be consistent with thermodynamics 
and quantum theory. We thus split the polarization and magnetization operators into 

P(x,t) ^ P fl (x,i)+P ext (x,i) 

M(x,t) ^ M fl (x,i) + M ext (x,i) (52) 

where the terms with the subscript 'fl' describe the fluctuations in the material and the 
'ext' all other sources like the dipole moments of atoms or molecules. In thermal equi- 
librium, the fluctuations average to zero, as we found after Eq. ll28l . In the following, 
more explicit information about the polarization noise is not needed. We shall assume 
that a Hamilton operator H exists that generates the macroscopic Maxwell equations 
as the Heisenberg equations of motion for the electric and magnetic field operators. 
(We are actually adopting a quantum Langevin picture, see Mandel and Wolf ] 12 [ and 
Gardiner 1 15 [.) A similar demonstration has been given by Wylie and Sipe 1 6 [. 
The field fluctuation spectrum, from Eq. is given by the expectation value 

+oo 

£(x,x';w) = J dre iwr (E(x,0)E(x',r)). (53) 

— oo 
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We assume thermal equilibrium without external sources and have used the station- 
arity of the correlation function. We now connect this spectrum to the linear response 
of the field to an external dipole oscillator, following Callen and Welton |9|. 

The solution for the electric field operator in the presence of a polarization source 
can be represented in terms of the Green tensor in the quantum theory as well because 
the substitution I l52t preserves the linearity of the macroscopic Maxwell equations. We 
thus get the field due to an operator-valued source, plus a term describing the free 
evolution of the field 

E(x,t) = Efree(x,t) 

+ |^e^|dVg(x,x / ;a;). [P fl (x';u,) + P ext (x';u,)] . (54) 
V 

There is a similar contribution from the magnetization that we do not need for the 
present discussion. In equilibrium, the free field operator has zero average, and we 
get the expectation value 

(E(x,t)> = J ^e- i ^|dVg(x,x / ;o;) • (P ext (x';a;)). (55) 
v 

The Green tensor can thus be identified with the linear response of the average field 
to a classical external polarization source (where P ext (x'; ui) is c-number valued). 

The linear response of the field can also be calculated directly from the Heisenberg 
equations. This provides us with an alternative expression for the Green tensor where 
equilibrium correlations will become apparent. For simplicity, we focus in the follow- 
ing on the response to an electric point dipole at the position x'. The coupling of the 
field to the dipole is described by adding to the Hamiltonian the term 

if int (t) = -d(t)-E(x',t), (56) 

and the Heisenberg equation reads 

^E(x,i) = -£[E(x,t), ff]+i[E(x,i), £,(x',t)]c^), (57) 

where summation over j is understood in the last term. The first term generates the 
free evolution of the field. Solving to first order in d and taking the average, we 
identify the field response function in the time domain as 

+oo 

(Ei(x,t)) = J drx^xV^-r) (58) 

— oo 

Xi , ( x,x',r) = ( im^),E^,t-r)]) forr>0, ^ 
[ for r < 0, 
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where the time dependence of the field operators is that of the evolution under H. 
The response function is thus itself a correlation function of the field. Due to the 
stationarity of equilibrium, Eq. (59i does not depend on t. It can be checked directly 
that Xr/ ( x ) x '; r ) is rea l as in the classical theory. Taking the Fourier transform of Xij/ 
we thus get an expression for the Green tensor in terms of a field correlation spectrum 

+00 

&i(x J x» = ~ |dre iwr ([^(x,r), £,(x',0)]>. (60) 


By causality the time integral is running over one half of the real axis only. We can 
make an integral over all r appear, as it occurs in the fluctuation spectrum ll53t . by 
forming the combination 

— {gji(-x',x;uj) - [^(x,x';u;)]*} = 

+00 

Jdre^dE^O),^'^)]). (61) 

—00 

The following relation allows to permute operators occurring in equilibrium correla- 
tion functions: 

+00 +00 

J dTe iu)T (A(T)B(0)} = / dr e iu)T (B{0)A(T)). (62) 

—00 —00 

In the classical theory, h = and operator ordering is irrelevant. Using this identity in 
the second term of the commutator in Eq. ll6TTl . we find the FD theorem: 

£, M '.^) = e ^ 1 gj ' (X ' Xi " ) ^ gij(X ' X>)1 ' - <03> 

The form J4H is recovered when the Green tensor satisfies the symmetry relation 1I2TI . 
This requires the additional assumption that permittivity and permeability are sym- 
metric, an assumption that we shall make in this contribution. 

Eq. can be proved using the Gibbs density operator i2~6l and the solution for 
the Heisenberg operator A{r) [see also Eq. i2"Hlll: 

Tr \exp(iHT/h)AeM-iHT/h)BeM-f3H) 

e^(A(r)B(0)) = e^—± L (64) 

Tr exp(— pH) 

One shifts the integration path in the complex r-plane to the line —00 — ihf3 . . . +00 — 
ih(3 and assumes that for |r| — > 00, the correlation function vanishes (otherwise this 
limiting value can be subtracted). Along the shifted path, Eq. i64l becomes 



Tr 



exp((3H) exp(iHT/h)Aexp(-iHT/h) exp(-(5H)B exp(-/?#) 



Tr exp(-0H) 
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e e 



Tr Bexp(WT/h)Aexp(-iHT/h)exp(-pH) 

JShuj iujt L 



Tr exp(-f3H) 

= e l3huJ e iu)T {B(0)A(T)), (65) 

using cyclic permutation under the trace. The r-integral now yields the right hand 
side of Eq. ll62t . 

1.3.4 Non-equilibrium situations 

A typical non-equilibrium situation that occurs in physics on the nanometer scale is a 
temperature gradient inside a nanostructure. In thermal scanning probe microscopy 
to quote another example, structures are held at different temperatures, being in con- 
tact with different reservoirs. These kind of settings can be described by a slight gener- 
alization of the present theory provided one assumes that each volume element of the 
solid structure is locally in thermal equilibrium at the temperature T(x). In this case, 
we can write down the fluctuation dissipation theorem for the thermal polarization 
field, by generalizing Eq. i47l , 

t>( ' \ 2ftllm£ £(x;w) . 

V{x,x;u) = — -«5(x-x). (66) 

exp[huj I k-qI (xJJ — 1 

We have assumed a local dielectric response for simplicity. The corresponding field 
fluctuation spectrum can be computed from the field operator <54t where the freely 
evolving field (with material damping, but without material fluctuations) and the 
Green tensor appear. Without external sources, one gets 



£ij(x,x';a;) = £^ ree) (x,x ;w) (67) 

2 film eoe(r;u) 



+ Jd 3 r [G ifc (x,r;w)]*G ife (x',r;^ 
v 



exp[huj/k B T(r)] - 1' 



See Henry and Kazarinov for a similar approach 1161 . The first term is nonzero for 
a bounded material surrounded by a non-absorbing dielectric (like free space) and 
describes the photons incident from infinity towards the observation points x, x'. It 
accounts for all of the field fluctuations when there is no material absorption at all. If 
the field in the surrounding medium as assumed to be at zero temperature (like in the 
visible frequency range), this term is zero for uj > 0. Even at finite temperature, how- 
ever, this term is typically negligible at sub-wavelength distances from an absorbing 
structure. Under these conditions, the second one dominates i.e., the radiation due to 
the polarization noise {66j|. We summarize explicit examples above planar substrates 
in the next Section. 
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Example: thermal van der Waals force. Consider a planar, dielectric substrate 
at temperature Ti/ks > 0, with the empty half-space above at zero temperature. This 
describes a thermal source surrounded by a vacuum chamber with absorbing walls 
at zero temperature |16|. (We neglect, of course, the heating of the walls due to the 
radiation from the substrate.) This model is the idealization of a situation where the 
dielectric is heated to a temperature higher than its surroundings. 

The fluctuations of the electromagnetic field now have two contributions: a first 
one coming from the thermal currents in the substrate, and a second one coming from 
the vacuum fluctuations in the empty half-space. Written schematically, 

(Ei(r;u;)Ej(r';uj')) = 2tt5(u - u') {Wy[w;T 1 ;r,r / ] +^>;T = 0;r,r']} . (68) 

There is no cross term because the two sources are not correlated. The spectrum Wij 
describes the radiation of the substrate and is computed using the Lifshitz model: it is 
the radiation produced by thermal current fluctuations inside the substrate. The spec- 
trum Vij describes the vacuum field in the empty half-space. It may be calculated a la 
Lifshitz by allowing for a nonzero imaginary part into the vacuum dielectric constant 
that is put to zero in the final result | i6l ll7l . Alternatively, one can perform an explicit 
field quantization in the half-space 1 18 1 and retain only the modes incident from the 
vacuum side. The result is the same in both cases, of course. 

The contributions in Il68t can be combined in the following way to recover the 
equilibrium situation at T\ = 0. We recall that for oj < 0, the substrate field spectrum 
is proportional to 1 + n(\oj\/T\), while the vacuum spectrum is proportional to 1 + 
= 1. Similarly, for positive frequencies, we have Wij(u)) oc n{uo/T\), and Vij = 0. 
Schematically, we may write 

Wij&Tr, r,r'] = W ij [u;T= 0;r,r'}+n(\uj\/T 1 )W ij [\u;\;T = 0;r,r'] (69) 

where the second term vanishes at zero temperature and is even in to. Those terms 
in ll68l that survive at T\ = combine to give the zero temperature fluctuation- 
dissipation theorem. 

The radiation emitted by the substrate (the second term in ll69l ) gives the explicitly 
temperature-dependent radiation. Note that this term is not the imaginary part of the 
Green tensor because we are not dealing with an equilibrium situation. Its behaviour 
as a function of frequency and the interpretation of the corresponding atom-surface 
force is given in 1 19 [. 

1.4 Example: planar surface 

In this section, we calculate the properties of the electromagnetic field close 
to a planar surface. The fluctuation-dissipation theorem is used, and a plane 
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wave expansion of the Green tensors for the electric and magnetic fields is 
worked out. Power law asymptotics are determined. This text is adapted from 
the habilitation thesis "Coherence theory of atomic de Broglie waves and elec- 
tromagnetic near fields" (C. Henkel, Universitat Potsdam, April 2004). 

1.4.1 Electric and magnetic Green tensors 

Consider a nonmagnetic solid with (relative) permittivity e and permeability fi = 1 
that fills the half-space z < 0. We shall be interested in the fluctuation spectrum of 
the electromagnetic field in the vacuum half-space, in particular for sub-wavelength 
distances < z < A = 2-jtc/uj. Up to hundreds of nanometers from the surface, this 
regime is relevant even at optical frequencies. 

For the Green tensor C/(x, x'; uS), we can make the ansatz, provided both x and x' 
are located outside the solid, 

0(x, x'; u) = (x, x'; u) + (x, x'; u) (70) 

where C/( vac ) is the vacuum Green tensor. C/( rcfl ) describes the electric field reflected 
from the solid and is determined from the boundary conditions for the electric and 
magnetic fields at z = 0. The decomposition l l70l is convenient to compute the field 
fluctuation spectra via the FD theorem l|4lTl because it exhibits clearly the additional 
contribution due to the scattering from the surface |6 . 5|. This statement remains true 
for scatterers of arbitrary shape, with more a complicated expression for the reflected 
or scattered field, of course. 

Above a planar solid, the so-called Weyl expansion or angular spectrum represen- 
tation provides a natural plane-wave basis for the incident and reflected fields, see the 
textbook by Nieto-Vesperinas |20|. The xy-plane naturally plays a distinguished role 
here. Introducing two-dimensional in-plane wave vectors Q = (q x , q y ), we shall use 
the notation 

q(±) = Q±mfe. (71) 

where n is the unit normal. One then has the following Fourier expansion for the 
reflected Green tensor |6| 

z,z'>0: g( rcfl )(x,x';cj) = (72) 
iiv 2 f d 2 Q e i q(+)' x - i q(-)- x/ _ 
^JWf q z ^ p r ^ t(+)eM(_) ' 

where the integral runs over all in-plane wave vectors Q. The wave vectors for down- 
ward and upward waves (both propagating and evanescent) are given by q(— ) and 
q(+), respectively. The disc |Q| < uj/c corresponds to propagating waves where 



^(iv/cf - Q 2 , lmq z >0 (73) 
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is real, while | Q | > lo/c describes evanescent waves that decay or increase exponen- 
tially with distance. Evanescent waves are required in the Green tensor to describe 
correctly the near field of a point dipole. For the reflected field, they provide the dom- 
inant contribution at sub wavelength distances from the solid. All elementary plane 
waves satisfy the vacuum dispersion relation (q(±)) 2 = (lo/c) 2 . Their polarization 
vectors are given by 



e p (±) = ^±Z_^ = ±^£_^1 ) ( 74) 



e s (±) = Q x n, 

q( ±) x e s _ ±Qq z - nQ 

lo/c lo/c 

where Q is the unit vector along Q. We have normalized the polarization vectors i74l 
such that e 2 (±) = 1; note that no complex conjugation is involved, although the vec- 
tors are complex in general. The conventional polarizations s and p are also called 
TE and TM in the literature (TE = electric field transverse to the plane of incidence 
spanned by n and Q). 

The reflection from the solid mixes downward with upward waves and is charac- 
terized by the Fresnel reflection coefficients I2"ll l22l 



^e(lo)(lo/c) 2 -Q 2 (75) 



q z + ^e(uj)(uj/cy 2 -Q 2 
s{lo)(lo / c) 2 - Q 2 - e(uj)q z 



(76) 



y/e(u)(u/c)*-Q* + e(u>)q z 

Note that the permittivity of the solid only enters via the Fresnel coefficients. As long 
as the planar symmetry is not broken, Eq. J72t can also be used above a multilayer 
medium |23 24 1. Note that sign conventions differ for the Fresnel coefficients and 
polarization vectors; only the product of and the polarization vectors appearing 
under the sum in Eq. i72l has an unambiguous meaning. 

We define the magnetic Green tensor by analogy to Eq. i42l as the magnetic induc- 
tion field radiated by a point magnetic moment 

Bi(x;u) =Ti i j(x,x. f ;uj)fj 1 j. (77) 

From the Maxwell equations we find that in terms of its electric counterpart, the mag- 
netic Green tensor in the vacuum above the solid is given by the double curl 

1 d d 

Wij(x,x';w) = —^eikiejmn^ — -£//„(x,x ;cu). (78) 
LO 2 OTfc ox' m 

For the free space Green tensor, this leads to an expression similar to the electric case. 
For the reflected field, as given by the Green tensor S72\ , we observe that the curl 
exchanges the polarization vectors i74l according to 

LO LO 

q x e s = -e p , qxe p = e s (79) 

c c 
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because q, e S/ and e p form an orthogonal Dreibein. Hence, up to a factor 1/c 2 , we ob- 
tain the magnetic field reflected from the solid by exchanging the reflection coefficients 
r s r p in Eq. J721 . 

1.4.2 Short distance expansions 

To illustrate the behaviour of the field at short distances, we review here asymptotic 
expansions for the electric and magnetic Green tensors. 

Electric field. As a first step, we show that in the near field, the reflected part i72l 
of the electric Green tensor takes a simple, electrostatic form. The integral over the 
wave vector Q involves the factor e l iz( z + z ') which provides a natural cutoff for large 
Q as soon as q z becomes imaginary with \q z \ > l/(z + z') 3> oj/c. Analyzing the 
integrand, we notice that it peaks around the cutoff value. We thus get the leading 
order asymptotics by using an expansion for Q, \ q z \ much larger than u/c under the 
integral. 

Let us first assume the more stringent condition Q S> \y/e\u)/c which corresponds 
to a distance much shorter than the medium wavelength, z -C \/\t/e\. The reflection 
coefficients and polarization vectors then behave like 

r s e s (+)e s (-) « (e-l)^(nxQ)(nxQ), (80) 

r p e p (+)e p (-) « £ ~ ) Q2 f (n - iQ)(n + iQ), (81) 

where n is the surface normal. The polarizations behave very differently, the p- 
polarized part dominating for large Q by a factor (Qc/uj) 4 . The reflection coefficient 
r p tends towards the electrostatic value r sta t = (e — + 1) and becomes indepen- 
dent of Q. This means that the reflection from the surface is nondispersive and can 
be modeled in terms of image theory ]21 1: the reflected field corresponds to the well- 
known field of an image dipole d = (—d x , —d y , d z )r stait located at the position X' — z'n 
below the surface. Performing the integrals over Q, we indeed find the short-distance 
asymptotics 

g<p*)(X' + R,z,X',*?;u) (82) 
_ e(w) -1{R 2 + z 2 )l - 3RR + {z 2 - 2i? 2 )nn + 3z (Rn - nR) 
~ e{uj) + 1 47re {R 2 + z 2 f' 2 ' 

where z = z + z' and (R 2 + z 2 ) 1 / 2 is the distance between the observation point and 
the image dipole. This expression depends on frequency only via the electrostatic 
reflection coefficient. Note also the broken isotropy of the correlation tensor where 
the coordinates parallel and perpendicular to the surface appear in non-equivalent 
ways. Nevertheless, the symmetry relation $T\\ is satisfied. 
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In a similar way, an asymptotic expression for distances larger than A/|y^i| can be 
worked out. This is particularly interesting above metallic surfaces where e « 'kj/equj 
can be very large. In this case, the expansion l|82t is valid for z -C 5^ where the skin 
depth 



'2e c 2 



(83) 



can be much smaller than the wavelength (typically, at frequencies below the in- 
frared). For the complementary range 5 U -C z, one finds, repeating the analysis lead- 
ing to the asymptotics ll80l , i8ll (see |25[) 



ImC/(X, z, X, z; u>) 



to 2 6^1 

327T£oC 2 Z 2 



(84) 



For simplicity we give only the imaginary part of the tensor with coinciding positions. 
Note the different power law with distance z and also the isotropic noise strength 
regarding the field polarization. 

To summarize, the electric field fluctuation spectrum derived from the Green ten- 
sor ll82l using the FD theorem d41~b is 



£ (x, x; oS) 



1 + nn 



1 



167re (e^ - 1) 



S w z 2 

16wl 

'^61 



for z < 5u, 

for 8^ < z <C (<LA) 1/2 , 
for (^A) 1 / 2 < z. 



(85) 



It is only at distances z > (S^X) 1 ^ 2 that the blackbody spectrum, originating from 
the free space contribution (third case in Eq.i85t). becomes the dominant contribu- 
tion. Closer to the surface, the electric field fluctuations have a noise spectrum that 
can exceed the Planck formula by several orders of magnitude. Note also that the 
low frequency limit is given by a constant spectrum oc k^T/ (a z 3 ). The electric field 
fluctuations thus behave like white noise on the nanometer scale. 



Magnetic field. Analogous calculations give for the Green tensor of the magnetic 
field at distances below the skin depth, z,z' <C S u , 

ft( refl )(X' + R,z,X.',z';u) (86) 

~ T^? ( f^j ~ \ [r 1 - RR/(r + z) + znn + (Rn - nR)] 
Aire 2 [4r(r + z) 

£(u) - 1 z 1 + RR/(f + z) - znn 

e(lo) + 1 r(f + z) 
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where f = (R 2 + z 2 ) 1 / 2 . From the fluctuation-dissipation theorem, we get the fluctua- 
tion spectrum of the magnetic field. We give here also the regimes of larger distances 

1 + nn 



o(x, x; u) 



for z <C 5 W , 

z 

81(1 + nn 



8?r5 2 (e^ - 1) 



3^4 

1 



for <L « z < (^A 3 ) 1 / 4 , (87) 

, for (^A 3 ) 1 / 4 < z. 

Note the different exponents for the power laws with distance compared to the electric 
field, and the larger cross-over distance to the blackbody radiation spectrum. The 
low-frequency limit of the magnetic noise spectrum is oc ^ksTa/z, it is frequency- 
independent as well. Similar expressions have been derived in I26l l2"5ll27l . 

1.4.3 Electromagnetic energy densities 

Some of the results summarized above have been discussed by Joulain, Carminati, 
Mulet, and Greffet in a recent paper on the definition and measurement of the LDOS 
close to planar surfaces |28[. These authors analyze the spectrum of the electric and 
magnetic energy densities u c (z;uj), u m (z;uj) as a function of distance and point out 
that the definition of the LDOS should be taken with care given the non-equivalent 
role played by the electric and magnetic fields in the near field. To illustrate this, we 
plot in FigureQ]the ratio u m (z;u)/u e (z; lo). Notice the strong dominance of the mag- 
netic energy throughout the near field range up to z ~ A for a metallic surface. The 
asymptotic formulas l|85l . ll87t provide good agreement with a numerical calculation 
based on the exact plane wave expansion \72\ for the electric Green functions and its 
magnetic equivalent. 

For a dielectric surface with an essentially real permittivity, the near field energy 
is dominantly electric, as shown in Fig. ^ (right). This behavior is not covered by the 
asymptotics {85j, ll87l because the assumption |e| » 1 breaks down, but can be found 
from ll82l and its magnetic counterpart. At distances comparable to the wavelength, 
oscillations appear for both metallic and dielectric surfaces, that correspond to the 
standing waves formed by the partial reflection of the field. In the far field regime 
z> A, the symmetry between electric and magnetic energy is restored, as expected 
from the Planck formula. 



1.5 Discussion of the van der Waals force 
1.5.1 Short distance limit 

The dynamic polarizability for an atomic ground state has been given in (|10)| . 
A transition frequency uo eg provides a natural scale for the cu-integral of the van 
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Figure 1: Ratio of magnetic to electric energy density vs. distance from a medium- filled 
half-space, normalized to the wavelength A. Left panel: metal with e = 1 + 800 i (skin depth 
5 = 0.008 A). Right panel: dielectric (e = 2.3 + 0.1 i, S > A). The dashed and dotted lines 
correspond to the asymptotic expansions 1851 , 1871 . 



der Waals force. As a first estimate, we can assume that frequencies around 
the strongest electric dipole transition, uj ~ uj eg contribute most strongly to the 
integral. This allows to define the short distance limit: the distance is small 
compared to the corresponding (reduced) wavelength, z <C X eg = c/uj eg . Note 
that the same physical limit is achieved by neglecting retardation, c — > oo, or 
focussing on low frequencies, uj ~ uj eg — > 0. 

Now, the low frequency (small distance) limit of the Green tensor is simply 
given by electrostatics. For a planar surface, the reflected field can be con- 
structed in terms of an image dipole located below the surface at a symmetric 
point. From this construction, one finds a Green tensor 



,( SC )/ . \ 1 - I5ij + 5 iz 5 jz 



4ire E(uj) + l (2z 



1 3 



(88) 



as also outlined around Eq. (j82)> . The van der Waals interaction energy is now 
given by 

Kdw(^) « -^4^E(l(^l d |e)| 2 +l^l4|e)P 



2ne {2z) 

x Im 7 ^) -l u 

J 2tt e(u) + 1 u* - uj 2 - Oiu; v ' 

This expression is used to interpret high-resolution atomic spectra near sur- 
faces, see the lecture by D. Bloch. Particularly interesting features occur when 
the atomic resonance, uj ~ uj eg , gets close to a "surface resonance", e(oj) ~ — 1. 
If a model for the dielectric function valid at complex frequencies is known, 
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one can also evaluate the integral in d89b along the imaginary axis: 

(90) 







27re(iO + l«& + £ 



All quantities here are real and positive so that no real or imaginary part is 
needed any more. 4 

1.5.2 Large distance limit 

At large distances, we have to use the full Green tensor. It can be written in the 
form 

<?H(r, r; w) = — J KdK — A« r s (e) + (k% z 5 jz - 7 2 A^) r p (e) 

(91) 

where A^- = diag(l, 1, 0) in a frame where the z-axis points normal to the sur- 
face. The distance occurs in the exponential 



exp(2i 7 z) = exp(-2z^K 2 - lu 2 /c 2 ). (92) 

One shifts the ^-integration onto the imaginary axis and observes that in the 
limit z — > oo, only small values of K and u = i^ contribute. We can thus 
expand the integrand for small K and £ and evaluate the integral. We find 
to lowest order the static atomic polarizability a^ tat \ The combination k = 



K 2 + £ 2 /c 2 is a good integration variable, and the integral gives the famous 
1/ z A power law of the "Casimir-Polder" interaction. The remaining integral 
can be written in terms of an angle a (with £ = k cos a), and we get 

(stat) 



3hca\ 
ln 2 £o(2zy 

tt/2 



^cp(^) = - _ 9 , l_ M x (93) 



x 



J da sin a{ - cos 2 a A {j r s (£ sta t) + (2 sin 2 a 5 iz 5 jz + Ay) r p (e sta t)}. 
o 



The reflection coefficients are evaluated with the static dielectric constant, but 
depend in general on a, for example: 



1 — V £stat cos a + sin a 

r s (£ sta t) = V , (94) 

1 + V £ s t a t cos 2 a + sin 2 a 



4 A response function at u> = i £ corresponds to the response to a real-valued excitation that 
is exponentially switched on, oc e^'. Since the physical response is a real (hermitean) variable, 
the response function is real as well. 
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For a perfect conductor and also for metals with a nonzero DC conductivity, 
e sta t = ioo, and therefore r s = — 1 and r p = +1. In that case, the a-integral 
gives §<%. 

2 Transitions 

In this second part of the lecture, we discuss transitions between different 
quantum states of a trapped atom that are induced by electromagnetic fluc- 
tuations. For atoms trapped close to surface, the relevant quantum states are: 
Zeeman or hyperfine sub levels that are trapped in potentials of different steep- 
ness (or not trapped at all), and center-of-mass eigenstates in the trapping po- 
tential. In the first case, transitions lead to trap loss or to the decoherence of 
Zeeman or hyperfine superposition states. In the second case, field fluctua- 
tions feed energy into the center-of-mass motion and give heating. 

The transition rates can be computed, as a first approximation, with 
second-order perturbation theory, using Fermi's Golden Rule. We shall find 
out that the rates are indeed small so that higher order corrections are usually 
not needed. We show in particular that the transition rate is a meausre of the 
cross spectral density of the field fluctuations. Let us look first at the transition 
between different magnetic sublevels. The following section is adapted from 
the paper "Fundamental limits for coherent manipulation on atom chips" by 
C. Henkel, P. Kriiger, R. Folman, and J. Schmiedmayer [Appl. Phys. B 76 (2003) 
173]. 

2.1 Noise spectrum and spin flip rate 

The coupling of the atomic magnetic moment to fluctuating magnetic fields gives rise 
to both spin flips and changes in the center-of-mass motion (scattering). The rate of 
these processes is given by Fermi's Golden Rule. We recall here that it can be con- 
veniently expressed in terms of the noise spectrum of the magnetic field fluctuations. 
(See |6[ for a similar approach and Chap. IV of ]29| for the derivation of a full master 
equation.) 

If we write | i) and | f ) for the atomic states before and after the transition, the tran- 
sition rate is 

IUf = |<F,f|tf int |I,i>| 2 <5(£ F + E t -E 1 - E,), (95) 

11 F,I 

where |I) and |F) are initial and final states for the field, the summation being an 
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average over the initial field states (with probabilities p(T)) and a trace over the final 
field states. The interaction Hamiltonian is given by H int = — fj, ■ B(x). 

Consider first the rate for spin flips. Since only a subset of magnetic sublevels 
\m\) are weak field seekers, spin flips \m\) — > |mf) are responsible for trap loss. The 
magnetic field is evaluated at the position r of the trap center. (An average over the 
atomic position distribution would be more accurate.) We write the ^-function for 
energy conservation as a time integral over e^i-^F-^fi)*/^ where Tiun = Ef — E\. 
The exponential e 1 ^ 1- - 5 ^*^ can be removed by introducing the field operators in the 
Heisenberg picture 

B(r, t) = e iH t/h B ^ e -iH t/h (%) 

(here, H is the free field Hamiltonian) and taking matrix elements of this operator 
between the initial and final field states. This gives 

e i(£l - Bp)t/n (I|B(r)|F) = (I|B(r,t)|F). (97) 

The sum over the final states |F) now reduces to a completeness relation and we get 
(a, (3 denote field components) 

27rn^p(I)(I|S Q (r)|F)(F|^(r)|I)«5(S F -Eh- hu) 

F,I 

oo 

= / dte^ t J2pm\B a (r,t)B f3 (r,0)\l) 

= B a p(r,r;u). (98) 

In the last line, we have defined the magnetic noise spectrum which is the Fourier 
transform of the field's autocorrelation function. The rate for spin flips can now be 
written as 

Ti^f = ^2 ( m i\f 1 a\m { }('m i \fJ'i3\m i }B a i3{r,r,-uj n ). (99) 

a,/3=x,y,z 

Since rrif / mi, the matrix elements of /i Q are only nonzero for directions perpen- 
dicular to the magnetic field at the trap center. We also recover the selection rule 
rrif — m\ = ±1 so that the relevant transition frequency is the Larmor frequency 
l^fil = <^L- The spin flip rate gives the order of magnitude of trap loss even if more 
than one weak-field seeking Zeeman states, m\ = +2, +1, say, are trapped (possible 
with many of the alkali atoms). This is because the matrix elements between adja- 
cent sublevels do not significantly differ in magnitude so that the atoms reach the 
non-trapped sublevel ni{ = after a time ~ 2/r + 2_>+i. 

We finally note that as long as the behaviour of the 'environment' (the field) is ig- 
nored in the description of the atom's dynamics, the noise spectrum is the only quan- 
tity needed to characterize the environment. It is also an experimentally measurable 
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quantity: for example, the rms magnetic noise (-B^(r)) 1 / 2 measured by a spectrum 
analyzer in a given frequency band Aui/2ir around to is (2B xx (r, r; u) Aiv /2n) l l 2 , the 
factor 2 accounting for the sum over positive and negative frequencies. The atomic 
spin flip rate may be regarded as an alternative way to measure the noise spectrum. 
In order of magnitude, the magnetic moment is comparable to the Bohr magneton, 
{hb/2ttTl = 1.4MHz/G), and we get 

IW(r) ~ O-OIs^Wmb) 2 ^^ - (100) 

Note that current SQUID magnetometers are able to detect magnetic field noise even 
on the WfT/V^z scale |27|. 



2.2 Magnetic near field noise 

Eq. ((99)) clearly shows that we need the spectrum of the magnetic field fluctua- 
tions for an estimate of the spin flip transition rate. The fluctuation-dissipation 
theorem provides the spectrum in terms of the Green function Wy(x, x'; uj) for 
the magnetic field, 

%(x,x';w) = /T _ i Im^, J (x,x». (101) 

The regime typical for trapped atom experiments is sub-wavelength distances, 
where the magnetic Green tensor (remember: the field radiated by a point 
magnetic moment) is dominated by the part reflected or scattered from the sur- 
face. From the short-distance expansion dSoT l, we find two regimes, depending 
on the relatize magnitude of the distance compared to the skin depth inside 
the surface material. 



Very short distance, normal conductor. The skin depth 5 U is large compared 
to the distance. As a scaling law for a normal metal, 

^ = ™^{z0iiiT- (W2) 

where g is the specific resistance. The scaling of the spin flip rate is then 

r^ f oc — , (103) 

QZ 

This estimate applies to a half-space. For absolute numbers, see the lecture 
by J. Schmiedmayer. Note in particular that the flip rate is independent of the 
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Larmor frequency. This is due to a conspiracy between the "density of states" 
of the magnetic near field (oc to) and the Bose-Einstein occupation number 

(oc 1/uj). 

The estimate l|103|) can be extended to other geometries using the follow- 
ing argument. In this "very short distance" regime, the distance z gives the 
depth of the material that contributes significantly to the magnetic noise out- 
side the surface. This has been checked by an explicit calculation in Il27ll30ll . 
Any structure thicker than z thus essentially behaves like a half-space. One 
gets a reduction of the magnetic noise spectrum for thinner films or wires, for 
example. Each dimension along which the metallic material becomes "finite" 
(characteristic length a like thickness or wire radius) gives an additional fac- 
tor a I z for the noise spectrum and hence for the spin flip rate. It is likely that 
this explains the cross-over seen in the lifetime data vs. distance of J. Schmied- 
mayer's group. 

If one works in the short-distance regime, the optimal material choice is an 
insulator (large resistance g). For strong magnetic fields, permanently magne- 
tized structures can be used. Alternatively, materials may exist whose conduc- 
tivity is reasonably large in the DC range and sufficiently small in the kHz to 
MHz range that applies for the relevant transitions. 

Intermediate distance, very good conductor. At larger distances or equiva- 
lently, smaller resistance, the skin depth ultimately becomes smaller than the 
distance. In this regime, a different power law is obtained. One gets a steeper 
scaling with distance, 

IW cx - * , (104) 

u L z 

which is now also frequency dependent. This regime should also apply to su- 
perconducting substrates. Note that in the relevant frequency range, even a 
superconductor shows some finite resistance due to the normal electron frac- 
tion not bound into Cooper pairs. 

In this regime, one can say that qualitatively, only the surface of the mate- 
rial produces magnetic field fluctuations. The power laws for different geome- 
tries therefore do not probably show the cross-over seen at short distance. The 
problem is currently under study. Estimates for magnetic noise above super- 
conducting substrates are given in the review by J. A. Sidles |26J. 
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2.3 Other rates 



The rate estimates for other processs like heating (excitation of higher trap 
eigenstates) and decoherence (decay of a superposition of magnetic sublevels 
into a statistical mixture) are discussed in the review paper by R. Folman 
and colleagues |31 [ and in the review by C. Henkel and the J. Schmiedmayer 
group II32I . The most robust choice for superposition states are magnetic sub- 
levels that have the same magnetic moment. They are shifted in the same way 
by magnetic fluctuations so that their relative phase is preserved. Such states 
exists in alkali atoms, they have the same orientation of the electron spin and 
opposite nuclear spins. The difference in magnetic moment is only provided 
by the nuclear magneton which is much smaller than the Bohr magneton, by 
roughly the mass ratio between the nucleon and the electron. 

Conclusion 

Atoms are a sensitive probe of electromagnetic fields close to surfaces. Their 
energy level shifts depend on the reflection or scattering of the field from the 
surface, integrated over the entire frequency range. At short distances, the 
range around atomic transition frequencies gives the dominant contribution. 
This can lead, for a "hot" surface, to interesting non-equilibrium forces. 

Transitions between different "internal" atomic states or the heating in a 
trap close to a surface gives directly access to the field's cross spectral density 
at the transition frequency, provided perturbation theory holds. In magnetic 
traps that confine only a subset of Zeeman sublevels, the magnetic spectrum is 
the dominant driving agent for transitions and leads to trap loss, heating, and 
decoherence. The corresponding rates show different characteristic scalings 
with distance and substrate geometry, depending on the size of the skin depth 
inside the material. 
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